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cutting  at  great  depth,  primary  combustors  for  high-speed  airbreathing  propulsion  systems, 
afterburners,  fuel/air  explosions,  and  spark-ignidon  internal  combustion  engines. 

Studies  of  premixed  turbulent  hvdrogen/air  jet  flames  involved  combustion  within 
the  nearly  homogeneous  potential  core  region  for  jet  Reynolds  numbers  of  7(X)0-4(KKI(). 
Measurements  included  the  following:  laser  light  sheet  imaging  for  flame  surface  statistics, 
conditional  laser  velocimetry  for  gas  turbulence  properties,  Rayleigh  scattering  for 
unreactedness,  and  Hash  schlieren  photography  for  flow  visualization.  Predictions  of  flame 
structure  were  undertaken  using  both  contemporary  turbulence  models  and  stochastic 
simulation.  The  stochastic  Simula  ion  was  developed  during  the  present  study  and  involved 
use  of  statistical  time  senes  techniques  to  simulate  unburned  gas  velocities  along  the  name 
surface  combined  with  a  two-dimensional  flame  advection  and  propagation  algorithm. 
Major  findings  of  this  phase  of  the  investigation  are  as  follows:  ( 1 )  effects  of  preferential 
diffusion  are  relevent  for  flames  at  high  Reynolds  numbers,  retarding  and  enhancing  the 
distortion  of  the  flame  surface  by  turbulence  for  stable  and  unstable  conditions, 
respectively;  (2)  local  turbulent  burning  velocity,  flame  brush  thickness  and  the  fractal 
dimension  of  the  flame  surface  all  increase  with  distance  from  the  flameholder.  with  larger 
rates  of  increase  at  larger  turbulence  intensities;  (3)  estimates  of  flame  properties  using 
contemporary  turbulence  models  were  only  fair  because  these  methods  cannot  account  for 
effects  of  preferential  diffusion,  distance  from  the  flameholder  and  finite  laminar  flame 
speeds;  and  (4)  the  stocnastic  simulation  duplicated  measured  trends  of  flame  surface 
{'properties  for  neutral  preferential  diffusion  conditions  (the  only  case  considered  but 
underestimated  effects  of  turbulence  (particularly  near  the  flame  up)  due  to  the  limitations  of 
a  two-dimensional  simulation. 

Work  on  freely  propagating  flames  was  limited  to  neutral  preferential  diffusion 
conditions.  Experiments  were  earned  out  for  spark  ignited  flames  in  isotropic  turbulence 
within  a  fan  stirred  reaction  vessel  with  turbulence  Reynolds  numbers  in  the  range  1 9(X)- 
4200.  Measurements  included  the  following:  laser  light  sheet  imaging  for  flame  surface 
statistics,  two-point  laser  velocimetry'  for  unbumed  gas  turbulence  properties,  and  flash 
schlieren  photography  for  flow  visualization.  Predictions  of  flame  structure  were  limited  to 
the  two-dimensional  stochastic  simulation.  Major  findings  of  this  phase  of  the  investigation 
are  as  follows:  (1)  local  turbulent  burning  velocity,  flame  brush  thickness  and  the  fractal 
dimension  of  the  flame  surface  all  increased  with  time  (distance)  from  the  point  of  ignition, 
with  larger  rates  of  increase  at  larger  turbulence  intensities-  (2)  evi;;:ir.g  correlations  of 
turbulent  burning  velocities  and  flame  surface  fractal  dimensions  are  incomplete  because 
they  don’t  account  for  effects  of  time  (distance)  from  the  point  of  ignition;  and  (3)  estimates 
of  flame  surface  properties  based  on  the  two-dimensional  simulation  were  reasonably 
good,  using  the  artifice  of  doubling  the  velocity  fluctuations  to  account  for  the  limitations  of 
a  two-dimensional  simulation. 

Definitive  evaluation  of  the  stochastic  simulation  technique  will  require  extension  to 
a  three-dimensional  algorithm.  Fortunately,  extension  to  three  dimensions  is  tractable,  even 
for  high  Reynolds  number  flames,  due  to  the  computational  efficiency  of  the  velocity 
simuladon. 
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D3  fractal  dimension  of  flame  surface 
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Rej  Reynolds  number  of  turbulence,  u0'  HucAo 

Sl  laminar  burning  velocity 

Sy  turbulent  burning  velocity 

ST.eff  effective  turbulent  burning  velocity 

u  stream  wise  velocity 

v  crosstream  velocity 

x  streamwise  distance 


Greek  Symbols 

a  thermal  diffusivity 

5l  laminar  flame  thickness 

e  ruler  length  for  fractal  property  measurement 


vi  1  1 


V 


kinematic  viscosity 
fuel-equivalence  ratio 


0 


mean  burner  exit  value 
burned  and  unbumed  gas  property 
burner  exit  condition 


Subscripts 

avg 

b,u 

o 


Superscripts 

mean  and  r.m.s.  fluctuating  property 


1.  INTRODUCTION 


This  investigation  considered  the  properties  of  turbulent  premixed  flames  both 
theoretically  and  experimentally.  Attention  was  limited  to  combustion  of  mixtures  of 
hydrogen,  oxygen  and  nitrogen  because  of  interest  in  using  premixed  hydrogen/oxygen 
flames  for  underwater  metal  cutting  processes  at  great  depth  where  high  pressures  cause 
normally  laminar  cutting  flames  to  become  turbulent.  However,  the  research  also  has 
application  to  a  variety  of  other  premixed  turbulent  combustion  process;  for  example: 
primary  combustors  for  high-speed  airbreathing  propulsion  systems,  afterburners,  ‘uel/air 
explosions  and  spark-ignition  internal  combustion  engines. 

This  is  a  final  report  of  accomplishments  during  the  investigation.  The  present 
report  is  brief,  additional  details  can  be  found  in  other  articles,  papers  and  theses  generated 
during  the  course  of  the  work,  see  Table  1  for  a  summary.  The  articles  by  Wu  et  ai.  ( 1990, 
1991a,  b)  and  Kwon  et  al.  (1991)  are  appended  to  provide  a  convenient  reference  Finally, 
Kwon  (1991)  and  Wu  (1991)  provide  detailed  descriptions  of  the  research  as  well  as 
complete  tabulations  of  data. 

The  investigation  was  divided  into  two  phases,  one  considering  turbulent  jet  flames 
because  they  approximate  the  configuration  of  flames  used  in  flame-cutting  processes,  the 
other  considering  freely-propagating  flames  in  isotopic  turbulence  because  this  fundamental 
configuration  simplifies  the  interpretation  of  results  needed  to  gain  a  better  understanding  of 
turbulence/flame  interaction.  Highlights  of  each  of  these  phases  of  the  research  are 
considered,  in  turn,  in  the  following  closing  with  a  summary  of  the  major  conclusions  of 
the  study. 


2.  TURBULENT  JET  FLAMES 


2.1  Introduction 


Turbulent  premixed  flames  are  an  important  fundamental  combustion  problem  with 
numerous  practical  applications;  therefore,  they  have  received  significant  attention  in  the 
past,  see  Wu  et  al.  (1990,  1991a, b)  and  references  therein  for  review  of  past  work. 
Nevertheless,  experimental  and  theoretical  progress  in  the  field  has  been  limited.  The  main 
difficulties  are  that  existing  measurements  either  involve  complex  flows  that  are  difficult  to 
interpret  or  compute  accurately,  or  they  involve  moderat’  Reynolds  numbers  where  the 
degree  of  development  of  turbulence  and  the  intrusion  of  buoyancy  aie  not  representative  of 
turbulent  flames.  The  objective  of  the  present  phase  of  the  investigation  was  to  help  address 
these  limitations  by  measuring  the  structure  of  turbulent  hydrogen/air  jet  flames  at  high 
Reynolds  numbers  typical  of  practical  applications.  Hydrogen/air  flames  were  studied  due 
to  their  relevance  to  metal  cutting  and  propulsion  applications.  Hydrogen/air  flames  also  are 
of  interest  because  they  place  unusual  stress  on  concepts  of  modeling  turbulent  flames, 
e.g.,  they  have  unusually  high  laminar  flame  speeds  and  are  subject  to  diffusive-thermal 
instabilities  from  the  preferential  diffusion  of  hydrogen.  Thus,  to  help  quantify  the 
significance  of  these  phenomena,  the  measurements  were  used  to  evaluate  predictions  of 
the  process.  All  work  was  limited  to  the  thin  flamelet  regime  defined  by  Bray  (1980), 
where  the  flame  thickness  is  smaller  than  the  smallest  (Kolmogorov)  scales  of  the 
turbulence,  because  this  reginie  is  most  relevant  to  practical  applications. 

Hydrodynamic  and  diffusive-thermal  (preferential  diffusion)  instabilities  are  the  two 
main  mechanisms  of  instability  of  premixed  flames  (Clavin,  1985).  However, 
hydrodynamic  instability  is  weak  and  has  only  been  observed  in  special  circumstances 
(Groff,  1982);  therefore,  it  is  generally  agreed  that  preferential  diffusion  is  responsible  for 
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most  observations  of  laminar  flame  instability.  Within  the  thin  flamelet  regime,  preferential 
diffus:on  instability  involves  the  interaction  between  diffusion  of  a  faster-diffusing  r  actant 
and  the  variation  of  laminar  burning  velocity  with  fuel-equivalence  ratio  (Manton  :t  al., 
1952).  Reactants  are  consumed  within  a  chin  reaction  -> r »•  so  that  the  faster-diffusing 
reactant  accumulates  near  bulges  in  the  flame  surface  e;  .g  toward  the  reactants  and  is 
depleted  from  oulges  extending  toward  the  products.  as,  when  the  laminar  burning 
velocity  increases/decreases  with  increasing  concentration  of  the  faster-diffusing  reactant, 
the  bulges  grow/decay  and  the  flame  is  unstable/stable.  Hydrogen  is  the  faster-dil fusing 
reactant  in  hydrogen/air  flames,  while  the  maximum  laminar  burning  velocity  is  reached  at  a 
fuel-equivalence  ratio  (0)  of  1.8;  therefore,  these  flames  are  stable/unstable  for  0  greater  or 
less  than  1.8. 

Effects  of  preferential-diffusion  instability  have  largely  been  studied  to  explain  the 
instabilities  of  laminar  flames,  however,  they  also  affect  turbulent  flames.  Clavan  and 
coworkers  (1985)  have  studied  tins  pnenomenor,  for  turbulent  flames  having  large  scale 
and  low  turbulence  intensity.  They  found  tha»  tu.oulent  distortion  of  the  flame  surface  was 
chaotically  enhanced  tor  unstable  conditions,  while  the  flame  acted  like  a  high-pass  filter, 
damping  low-frequencv  deflections,  for  stable  conditions.  These  studies  involved  only 
weak  turbulence  however,  so  that  investigate  -  of  corresponding  effects  for  practical 
turbulent  flames  is  of  significant  interest. 

Exact  numerical  simulation  ot  p^mixed  turbulent  flames  is  too  computationally 
intensive  to  be  feasible  for  the  fo'^see-ble  future  (Pope.  1990);  'herefore,  turbulence 
models  and  approximate  numerical  simulations  have  received  significant  attention  (Bray, 
1980;  Pope,  1990).  Among  the  turbulence  models.  Cant  and  Bray  (1988)  report  a  relatively 
advanced  method  that  can  allow  for  flame  quenching  by  large  levels  of  turbulent  distortion 
'flame  stretch);  however,  this  approach  has  not  been  evaluated  using  measurements. 
Similarly,  approximate  numerical  simulations  have  been  described  that  models,  remove 
many  of  the  ad  hoc  features  of  turbulence  (Ashurst  and  Barr,  1983;  Ghoniem  et  al.,  1980, 
1983);  however,  these  methods  are  still  too  computationally  intensive  to  allow  their 
evaluation  using  measurements  from  large  Reynolds  number  turbulent  flames. 

The  objectives  of  the  ^resent  investigation  were  designed  to  help  fill  some  of  the 
existing  gaps  in  the  literature.  Measurements  were  undertaken  to  provide  the  turbulence 
properties  of  the  unreacted  gas,  and  flame  surface  statistics,  for  high  Reynolds  number 
turbulent  premixed  round  jet  flames.  These  results  were  used  to  assess  the  relevance  of 
preferential-d’ffusion  phenomena  for  high  Reynolds  number  flames  and  to  evaluate 
predictions  of  the  process.  Two  predictive  methods  were  examined:  (1)  the  contemporary 
turbulence  model  of  Cant  and  Bray  (1988),  and  (2)  an  approximate  numerical  simulation 
developed  during  the  present  investigation  that  has  sufficient  computational  efficiency  to  be 
tractable  for  high  Reynolds  number  turbulent  flarr.es. 

2.2  Experimental  Methods 

Apparatus.  The  test  apparatus  consisted  of  a  round  coaxial  burner  ’.here  the  outer 
burner  flame  provided  a  hot  environment  to  stabilize  the  inner  burner  flame  (which  was  the 
flame  that  was  studied).  The  burners  injected  vertically  upward  with  the  inner  burner  flow 
being  fully-developed  turbulent  pipe  flow  from  a  tube  having  an  1 1  mm  inside  diameter. 
The  burners  operated  at  atmospheric  pressure  with  the  exhaust  products  removed  through  a 
hood.  Instruments  were  mounted  rigidly  so  that  the  burner  was  traversed  horizontally  and 
vertically  to  access  various  points  within  the  inner  burner  flame. 

Instrumentation. Measurements  involved  flash  schlieren  photography  for  flow 
visualization,  laser  velocimetry  to  measure  gas  velocities  ana  turbulence  properties. 


Rayleigh  scattering  to  measure  the  extent  of  reaction  and  flame  tomography  to  measure 
flame  surface  statistics.  Flash  schlieren  photography  is  a  well-known  technique  so  only  the 
other  methods  will  be  discussed  in  the  following. 


A  dual-beam  forward  scatter  laser-velocimeter  (LV)  was  used  to  measure 
conventional  time-averaged  mean  and  fluctuating  velocities  in  the  flow,  as  well  as 
conditional  velocities  and  turbulence  properties  in  the  unreacted  gas.  Frequency  shifting 
was  used,  when  necessary,  to  control  directional  bias  and  ambiguity.  The  optical  plane  of 
the  laser  beams  was  rotated  to  measure  streamwise  and  radial  velocities.  For  overall  flow 
measurements,  the  flow  was  seeded  with  aluminum  oxide  particles  having  a  nominal 
diameter  of  1  pm.  Conditional  velocity  properties  were  found  by  seeding  the  flow  with  oil 
droplets  having  a  nominal  diameter  of  1  pm:  these  drops  vaporized  at  the  flame  front  so  that 
velocities  were  only  measured  in  the  unreacted  gas. 

The  Rayleigh  scattering  arrangement  was  similar  to  Driscoll  and  Gulati  (1988).  The 
488  nm  line  of  an  argon-ion  laser  was  focussed  at  a  point  in  the  flame  and  the  scattered 
light  from  the  point  was  observed  using  a  photodetector.  The  Rayleigh  scattered  light  from 
the  dense  unreacted  gas  was  much  stronger  than  from  the  combustion  products.  Thus,  the 
output  was  a  telegraph  signal  yielding  the  fraction  of  time  when  reactants  were  present  (the 
unreactedr.ess)  and  the  time-averaged  frequency  of  shifting  between  reacted  and  unreacted 
states.  All  other  statistical  properties  of  a  telegraph  signal  can  be  found  from  these  two 
properties. 

Flame  tomography  was  used  to  find  flame  surface  statistics.  This  involved  seeding 
the  reactants  with  oil  drops  having  a  nominal  diameter  of  1  pm  which  burned  out  in  the 
flame  front  similar  to  the  conditional  LV  measurements.  The  drops  were  illuminated  with  a 
2(X)  pm  thick  laser  light  sheet  from  a  pulsed  dye  laser.  The  laser  light  scattered  from  the 
drops  was  photographed  normal  to  the  light  sheet  with  the  short  duration  of  the  pulsed  laser 
dps)  acting  to  stop  the  flow  so  that  the  edge  of  the  drop  containing  region  (the  flame  sheet) 
could  be  resolved.  These  records  were  processed  using  an  image  analyzer  system  to  yield 
mean  and  fluctuating  flame  radius,  the  perimeter  of  the  flame  and  the  fractal  dimensions  of 
the  flame  surface  (see  Mantzaras  et  al.  (1989)  and  Mandelbrodt  (1982)  for  a  discussion  of 
fractal  properties  of  flame  surfaces).  Orienting  the  optical  system  in  the  vertical  and 
horizontal  directions  yielded  corresponding  flame  surface  properties. 

Test  Conditions.  Two  series  of  tests  were  conducted  yielding  a  wide  range  of  test 
conditions.  Table  2  is  a  summary  of  test  conditions  for  the  second  senes  which  was  used  to 
evaluate  effects  of  preferential  diffusion  and  the  performance  of  the  approximate  numerical 
simulation.  All  tests  were  in  the  thin  flamelet  regime  where  the  smallest  turbulent  length 
scales  (the  Kolmogorov  scales)  are  comparable  or  greater  than  the  flame  thickness.  Die  tests 
involved  4>  =  0.8,  1.8  and  3.6,  which  are  in  the  unstable,  neutral  and  stable  preferentral- 
diffusion  regimes.  The  unstable  and  stable  flame  conditions  involve  identical  laminar  flame 
speeds  and  jet  exit  conditions  so  that  differences  between  the  properties  of  the  flames 
provides  direct  evidence  of  effects  of  preferential  diffusion.  That  such  effects  were  observed 
can  be  seen  directly  from  the  flame  lengths,  Lc,  summarized  in  Table  2;  clearly,  L<-  is 
consistently  shorter  for  the  unstable  than  the  stable  flame  showing  that  stable  conditions 
reduce  distortion  of  the  flame  surface  by  turbulence  (and  correspondingly  the  turbulent 
flame  speed)  so  that  longer  flames  are  needed  to  complete  the  reaction. 


Table  2.  Summary  of  Turbulent  Jet  Flame  Tests3 


uo,avg/SL 

0.9 

2.5 

4> 

0.8 

1.8 

3.6 

0.8 

1.8 

3.6 

Preferential  Diffusion 

Unstable 

Neutral 

Stable 

Unstable 

Neutral 

Stable 

Regime 
uo,avg  (m/s) 

20.4 

31.1 

20.4 

55.0 

88.0 

58.3 

uo,avg  (m/s) 

2.0 

3.1 

2.0 

5.5 

8.8 

5.8 

Sl  (m/s)b 

2.2 

3.5 

2.3 

2.2 

3.5 

2.3 

Red=u0>avg<lA'o 

11300 

13900 

7000 

30600 

39200 

20000 

ReT=u  oavg  -^uc/Vo 

313 

389 

194 

861 

1 104 

562 

l</d 

2.6 

2.6 

2.8 

4.0 

5.3 

7.7 

flK(Bm) 

42 

35 

60 

20 

16 

27 

5L=oco/SL(jdm) 

19 

17 

35 

19 

17 

35 

Flame  regimec 

Wrinkled  thin  flamelets 

Mixing-limited  thin  flamelets 

aRound  jet  burner  (d  =  11  mm)  directed  vertically  upward  with  fiuo  =  3. 1  mm 
b Laminar  flame  speeds  from  Andrews  and  Bradley  (1973). 
cBased  on  definitions  from  Bray  (1980). 


2.3  Theoretical  Methods 

Turbulence  Model.  Theoretical  methods  evaluated  using  the  new  measurements  were 
based  on  procedures  of  turbulence  models  and  approximate  numerical  simulation.  The 
turbulence  model  involved  a  Favre-averaged  formulation  along  the  lines  of  Cant  and  Bray 
(1988)  with  some  features  drawn  from  earlier  work  in  this  laboratory,  see  Wu  et  al.  (1990) 
for  specific  details.  Major  assumptions  were  as  follows:  boundary- layer  approximations 
apply,  steady  (in  the  mean)  axisymmetric  flow  with  no  swirl,  low  Mach  number  flow  with 
negligible  potential  and  kinetic  energy  changes,  negligible  radiant  energy  exchange,  equal 
exchange  coefficients  of  all  species  and  heat,  high  Reynolds  number  flow  so  that  laminar 
transport  is  negligible  in  comparison  to  turbulent  transport,  negligible  effects  of  buoyancy 
on  turbulence  properties  and  infinitely-thin  flame  sheet  separating  the  burned  and  unbumed 
gases.  These  assumptions  are  either  representative  of  the  conditions  of  the  experiments  or 
are  typical  of  contemporary  turbulence  models  of  premixed  flames. 


Initial  conditions  for  the  computations,  at  the  burner  exit,  were  taken  to  be  fully- 
developed  turbulent  pipe  flow  having  properties  drawn  from  references  cited  in  Hinre 
(1975).  The  formulation  of  Cant  and  Bray  (1988)  accounts  for  quenching  due  to  turbulent 


distortion  of  the  flame  surface,  parameters  needed  for  this  approach  involved  laminar  flame 
speeds  from  Andrews  and  Bradley  (1973),  specific  recommendations  from  Cam  and  Bray 
(1988),  and  a  fit  of  present  measurements  at  one  test  condition  —  see  Wu  et  al.  (1990). 


Approximate  Numerical  Simulation.  Turbulent  flame  propagation  was  numerically 
simulated  using  statistical  time  series  simulation  of  the  velocity  field  of  the  unbumed  gas 
combined  with  a  flame  advection  and  propagation  algorithm.  The  simulation  of  the  velocity 
field  was  based  on  conventional  statistical  time-series  methods  (Box  and  Jenkins,  1976)  but 
extended  similar  to  earlier  work  in  this  laboratory  (Kounalakis  et  al.,  1991;  Parthesarathv 
and  Faeth,  1990)  to  treat  a  multidimensional  process,  see  Kwon  (1991)  for  a  derivation  of 
the  procedure  for  two-and  three-dimensional  time-dependent  velocity  fields.  The  flame 
advection  and  propagation  algorithm  was  taken  from  the  MIMOC  computer  program  of 
Ghoniem  et  al.  (1983),  which  is  limited  to  a  two-dimensional  time-dependent  simulation. 
This  is  a  major  limitation  for  treating  fundamentally  three-dimensional  time-dependent 
turbulence,  however,  assessing  the  present  approach  is  justified  before  undertaking  the 
significant  extension  required  to  treat  three-dimensional  effects. 

Other  major  assumptions  of  the  approximate  numerical  simulations  are  as  follows: 
constant-pressure  turbulent  deflagration  wave,  infinitely-thin  flame  sheet,  turbulence 
properties  of  the  unbumed  gas  are  stationary  and  are  unaffected  by  the  presence  of  the 
flame,  the  flame  is  neutrally-stable  with  negligible  effects  of  quenching  (i.e.,  relative  to  the 
unbumed  gas,  the  flame  propagates  normal  to  its  surface  at  the  laminar  burning  velocity). 
The  assumptions  of  constant  pressure,  thin  flames,  neutral  stability,  attachment  points  near 
the  burner  rim  and  turbulence  properties  unaffected  by  the  flame  were  all  conditions  of  the 
experiments  used  to  evaluate  these  simulations.  Neglecting  quenching  was  also  reasonable 
because  the  flames  were  at  the  maximum  laminar  burning  velocity  condition,  well  away 
from  flammability  limits,  and  values  of  the  relative  turbulence  intensity,  u’/Sl,  were 
relatively  small. 

Statistical  simulations  of  the  velocity  field  can  be  designed  to  satisfy  any  number  of 
the  statistical  properties  of  turbulence  but  priorities  must  be  set  to  control  computation  and 
computer  memory  requirements.  Thus,  present  calculations  were  limited  to  simulating  mean 
velocities,  Gaussian  probability  density  functions  of  velocity  fluctuations,  and  temporal  and 
spatial  correlations.  The  correlations  were  approximated  by  exponential  functions  because 
this  yields  a  Markov-like  simulation  that  vastly  simplifies  computations  (Box  and  Jenkins, 
1976;  Kwon  et  al.,  1991).  It  was  also  assumed  that  correlations  between  the  important 
nearest  neighbors  in  space  and  time  were  locally  homogeneous  and  could  be  decomposed  as 
products  of  correlations  in  each  coordinate  direction.  Under  these  approximations,  the 
simulation  of  velocity  fluctuations  at  a  particular  point  only  involves  results  at  the  seven 
nearest  neighbors  while  still  satisfying  correlations  over  all  space  and  time  (exactly  for 
homogeneous  isotropic  turbulence,  and  within  30  percent  for  locally-homogeneous 
turbulence). 

2.4  Results  and  Discussion 

Experimental  Observations.  Measurements  of  the  velocity  field  within  the  unreacted 
gas  indicated  that  flow  properties  approximated  turbulent  pipe  flow  at  corresponding  radial 
positions.  This  behavior  follows  because  the  present  flames  were  relatively  short  so  that  the 
unreacted  gas  was  largely  within  the  potential  core-like  region  of  the  flow.  Additionally, 
variations  of  turbulence  properties  were  not  large  over  this  region  and  approximated  locally 
homogeneous  flow  (Wu  et  al.,  1991a). 

Initial  flow  visualizations  using  flash  schlieren  photography  exhibited  effects  of 
preferential -diffusion  (Wu  et  al.,  1990).  Aside  from  the  flame  length  observations  discussed 


in  connection  with  Table  2,  the  flame  surfaces  for  stable  conditions  were  much  smoother 
than  for  unstable  conditions  when  other  properties  were  comparable.  This  implies  that  the 
preferential  diffusion  mechanism  acts  to  suppress/enhance  turbulent  distortion  of  the  flame 
surface  even  for  high-Reynolds  flames.  These  observations  could  be  quantified  by 
computing  average  turbulent  burning  velocities  for  the  entire  flame  surface,  Sjxff.  The 
resulting  values  of  Sj,  eff/^L  are  plotted  as  a  function  of  effective  relative  turbulence 
intensities  in  Fig.  1.  Values  of  Sl  in  these  correlations  are  drawn  from  the  measurements  of 
Andrews  and  Bradley  (1973).  Results  plotted  in  Fig.  1  are  for  0  =  0.8,  1.0,  1.8  and  3.6  for 
Red  in  the  range  7000-40000.  The  results  indicate  that  effective  turbulent  flame  speeds  fall 
on  two  branches:  a  neutral-unstable  branch  involving  0  =  0.8,  1.0  and  1.8  and  a  stable 
branch  involving  0  =  3.57.  The  stable  branch  has  substantially  less  response  to  increased 
turbulence  levels  (i.e.,  a  smaller  slope)  than  the  neutral-unstable  branch.  Furthermore,  even 
within  the  neutral -unstable  branch,  unstable  conditions  are  somewhat  more  responsive  to 
increased  turbulence  levels  than  the  neutral  condition.  These  trends  clearly  imply  that 
diffusive-thermal  phenomena  retard  turbulent  distortion  of  the  flame  surface  for  stable 
conditions  (thus  reducing  Sj.cff)  and.  to  a  lesser  degree,  enhance  distortion  for  unstable 
conditions  (yielding  higher  values  of  Sjcff)  for  high  Reynolds  number  turbulent  flames. 
After  finding  this  behavior,  review  of  earlier  measurements  in  the  literature  revealed  similar 
trends  for  other  fuels  like  H2  -  CH4  /air,  propane/air  and  iso-octane/air  flames,  as  discussed 
by  Wu  et  al.  (1991a),  see  Abdel-Gayed  et  al,  (1984)  for  a  typical  example.  Thus,  this 
behavior  is  not  solely  caused  by  the  unusually  high  mass  diffusivities  of  hydrogen  but  is 
also  characteristic  of  any  reactant  mixture  where  the  mass  diffusivities  of  all  species  are  not 
the  same. 

The  effect  of  both  preferential  diffusion  and  the  development  of  distortion  of  the 
flame  surface  with  increasing  distance  from  the  flameholder  can  bee  seen  from  the  results 
illustrated  in  Fig.  2.  This  involves  plots  of  flame  surface  area  ratios,  Aj/Al,  for  unstable, 
neutral  and  stable  conditions  (burner  exit  conditions  and  laminar  flame  speeds  are  identical 
for  the  unstable  and  stable  conditions,  0  =  0.8  and  3.6).  As  discussed  by  Gouldin  (1987), 
At/Al  is  proportional  to  St/Sl  at  at  modest  stretch  rates  for  neutral  conditions,  and 
approximately  so  for  unstable  and  stable  conditions.  First  of  all,  At/Al  is  seen  to 
progressively  increase  from  values  of  unity  at  the  flameholder  (corresponding  to  laminar 
combustion)  with  increasing  distance  from  the  burner  exit.  This  behavior  is  somewhat 
analogous  to  turbulent  dispersion  of  panicles  where  increased  distance  from  the  point  of 
introduction  of  the  particles  allows  the  action  of  turbulence  to  increase  the  lateral  spread  of 
the  particles.  Also  similar  to  effects  of  turbulent  dispersion,  the  rate  of  increase  of  Aj/Al  is 
increased  as  the  turbulence  intensity  increases.  Effects  of  preferential  diffusion  can  be  seen 
by  comparing  results  for  0  =  0.8  and  3.6  which  have  identical  values  of  u'/Sl-  The 
unstable  flame  generally  has  larger  values  of  Aj/Al  at  a  given  distance  from  the  burner 
exit,  reflecting  the  enhanced  distortion  of  the  flame  surface  for  unstable  conditions. 
However,  the  effect  of  growth  of  flame  surface  area  with  distance  from  the  flameholder 
plays  a  role  as  well;  for  example,  the  stable  flame  extends  farthest  from  the  burner  exit, 
due  to  lower  levels  of  distortion  near  the  burner,  but  eventually  exhibits  higher  values  of 
Aj/Al  than  the  rest  at  large  distances  from  the  burner  exit. 

A  popular  measure  of  the  degree  of  wrinkeldness  of  flame  tomographs  is  the  fractal 
dimension,  D2  (Gouldin,  1987)  which  can  be  measured  as  described  by  Mandelbrot 
(1982).  Results  obtained  from  horizontal  light  sheet  images  are  plotted  in  Fig.  3  for  the 
same  flame  conditions  as  Fig.  2.  Values  of  D2  are  seen  to  be  influenced  by  height  above  the 
burner,  u  0,avg/SL  and  effects  of  preferential  diffusion  instability,  similar  to  Aj/Al  in  Fig. 
2.  For  all  conditions,  D2  progressively  increases  with  increasing  distance  from  the  burner 
exit.  Near  the  burner  exit,  D2  approaches  unity,  which  is  the  value  representative  of  a 
smooth  grometrical  surface  (Mandelbrot,  1982)  —  corresponding  to  laminar  like  conditions 
mentioned  earlier.  This  follows  since  the  point  of  attachment  is  fixed  along  the  periphery  of 
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Fig.  1  Effective  turbulent  burning  velocity  of  jet  flame? 
as  a  function  of  relative  velocity  fluctuations. 
From  Wu  et  a  1 .  ( 1990 ) . 
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Fig.  2  Ratio  of  the  turbulent  to  mean  flame  area  as 
a  function  of  distance  from  the  flame  holder 
for  jet  flames.  From  Wu  et  al.  (1991a). 
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Fig.  3  Variation  of  fractal  dimension,  D^,  of  flame 

surface  with  distance  from  the  flame  holder  fo 
jet  flames.  From  Wu  et  al.  (1991a). 


the  burner  passage  exit  and  the  turbulence  does  not  have  sufficient  distance  (time)  to  distort 
the  flame  surface  from  this  boundary  condition.  The  subsequent  increase  of  D2  with 
distance  is  limited  by  two  factors:  the  finite  extent  of  the  flame,  so  that  the  largest  values  of 
D2  are  associated  with  the  longest  flames;  and  the  maximum  distortion  of  surfaces  in 
turbulent  flows,  which  appears  to  be  in  the  range  D2  =  1.3  -  1.4  (Gouldin,  1987; 
Mantzaras  et  al.,  1989,  Wu  et  al.,  1991a).  Notably,  results  for  the  longest  flame  in  Fig.  3 
appear  to  be  approaching  the  limiting  values  of  D2  in  spite  of  a  relatively  low  value  of 
u'o.avg/SL  and  the  fact  that  this  flame  is  stable.  Effects  of  preferential  diffusion  also  are  are 
evident  from  the  results  plotted  in  Fig.  3.  This  is  seen  most  readily  for  u'0>avg/SL  =  0.9 
where  the  initial  rate  of  growth  of  D2  is  more  rapid  for  $  =  3.6  (stable). 

Taken  together,  present  ooservations  of  the  development  of  flame  surface  distortion 
and  preferential  diffusion  preclude  notions  that  flame  surface  properties  on  based  solely  on 
local  conditions  and  only  involve  passive  propagation  of  the  flame  at  a  fixed  laminar  flame 
speed  (except  for  condidons  where  effects  of  quenching  are  small  at  the  singular  neutral 
preferential  diffusion  state).  Thus,  early  correlations  of  Sj/Sl  and  D2  (or  D3  =  D2  +  1) 
solely  as  a  function  of  uVSl  were  found  (Abdel-Gayed  et  al.,  1984;  Mantzaras  et  al.  1989) 
because  u/Sl  influences  the  rate  of  growth  of  these  properties  for  limited  ranges  of  flame 
sizes,  however,  effects  of  development  of  flame  surface  distortion  with  distance  clearly  are 
important  as  well.  Additionally,  preferential  diffusion  affects  of  the  rate  of  growth  of  these 
properties,  contribuung  to  apparent  scatter  of  such  correlations  even  for  fixed  values  of 
uVSl  and  flame  sizes.  Additionally,  the  growth  of  surface  distortion  with  distance  from  a 
flameholder,  even  when  turbulence  properties  remain  nearly  constant,  obviously  invalidates 
concepts  that  local  properties  alone  affect  turbulent  combustion  rates  —  ideas  which  are 
implicit  in  turbulence  models  of  combustion  that  have  been  proposed  thus  far.  That  these 
effects  have  been  observed  in  high  Reynolds  number  flames  implies  that  they  are  not  just 
artifacts  of  laboratory  experiments  but  must  be  considered  in  any  useful  theory  or  model  of 
turbulent  premixed  combustion. 

Theoretical  Results.  In  view  of  the  observations  that  were  just  discussed,  it  is  not 
surprising  that  the  comparison  between  predictions  using  the  Cant  and  Bray  (1988) 
turbulence  model  and  present  measurements  were  not  very  satisfying  (Wu  et  al.,  1990). 
Even  after  choosing  one  empirical  parameter  to  fit  measurements  and  predictions  at  a 
particular  condition,  the  model  was  not  very  effective  for  predicting  the  trends  of  the 
results.  The  intrinsic  difficulties  are  that  the  model  does  not  account  for  effects  of  distance 
from  the  flameholder  and  preferential  diffusion.  Another  problem  is  that  existing  turbulence 
models  implicitly  ignore  effects  of  laminar  flame  speeds  (implying  that  u'/Sl  very  large) 
which  is  not  appropriate  for  both  present  test  conditions  and  many  practical  applications. 
Finally,  although  the  Cant  and  Bray  (1988)  approach  allows  for  effects  of  local  quenching, 
these  effects  were  not  large  for  present  test  conditions;  therefore,  this  aspect  of  their 
approach  was  not  really  tested.  Taken  together,  contemporary  turbulence  models  will 
require  extensive  additional  development  to  become  effective  for  practical  turbulent  flames, 
in  spite  of  extensive  work  on  them  in  the  past. 

Due  to  the  limited  promise  of  turbulence  models  of  premixed  turbulent  flames,  the 
bulk  of  present  theoretical  work  was  devoted  to  approximate  numerical  simulations  (Wu  et 
al.,  1991b).  As  a  first  step,  considerations  were  limited  to  neutral-stability  conditions  and  a 
two-dimensional  time  dependent  simulation,  as  noted  earlier.  Simulated  flame  properties 
duplicated  measured  trends  of  the  variation  of  flame  surface  properties  and  mean 
unreactedness  with  distance  from  the  burner  exit  and  relative  turbulence  intensity. 
However,  effects  of  turbulence  were  underestimated  —  particularly  near  the  flame  tip  — 
due  to  the  limitations  of  a  two-dimensional  simulation  (i.e.,  due  to  out  of  plane  deflection 
of  the  thin  region  of  unreacted  gas  near  the  flame  tip).  The  simulations  of  freely 
propagating  flames  were  less  subject  to  this  problem  due  to  the  relatively  large-diameter 


flame  balls  that  were  considered;  therefore,  evaluation  of  the  approximate  numerical 
simulation  will  be  discussed  in  connection  with  these  results  in  the  next  section. 
Nevertheless,  present  computations  have  established  the  computational  feasibility  of  three- 
dimensional  time  dependent  approximate  simulations  —  even  for  high  Reynolds  number 
flames  —  so  that  the  method  clearly  merits  additional  development. 

3.  TURBULENT  FREELY-PROPAGATING  FLAMES 

3.1  Introduction 

An  experimental  and  theoretical  study  of  free  turbulent  premixed  flames  at  neutrally- 
stable  preferential  diffusion  conditions  is  discribed.  This  problem  is  of  interest  because  the 
constraints  of  flame  holding  and  the  complications  of  effects  of  preferential  diffusion  are 
absent,  while  the  flame  is  subjected  to  the  simplest  hydrodynamic  state  of  turbulence.  The 
experiments  involved  hydrogen/air/nitrogen  mixtures  ignited  within  a  fan-stirred 
combustion  chamber.  Measurements  included  flame  tomography  to  find  flame  surface 
statistics  and  two-point  laser  velocimetry  to  find  the  turbulence  properties  of  the  unbumed 
gas.  Test  conditions  were  in  the  wrinkled  thin  flamelet  regime  with  turbulence  Reynolds 
numbers  in  the  range  0-4195,  which  are  representative  of  the  intense  flames  found  in 
practical  applications.  Flame  surface  properties  were  numerically  simulated  using  time- 
series  simulation  techniques  (Box  and  Jenkins,  1976)  with  the  simulations  being  evaluated 
using  the  new  measurements. 

Earlier  measurements  of  premixed  flames  in  isotropic  tuibuience  have  been 
reported,  see  Abdel-Gayed  et  al.  (1984)  and  Kwon  ct  ai.  (1991)  for  a  discussion  recent 
work  along  these  line.  The  main  distinction  between  the  present  measurements  and  earlier 
work  is  a  greater  emphasis  on  the  evolution  of  measured  flame  surface  statistical  properties 
during  propagation  from  the  point  of  ignition  and  the  absence  of  preferential  diffusion 
effects.  Additionally,  it  is  hoped  that  the  small  rates  of  flame  stretch  in  comparison  to 
extinction  conditions,  the  simple  and  well-characterized  turbulence  field,  and  the 
information  on  flame  surface  development,  will  be  useful  for  developing  and  evaluating 
models  of  the  process.  As  discussed  in  connection  with  the  turbulent  jet  flames,  the  present 
approximate  numerical  simulation  represents  a  significant  departure  from  methods  reported 
in  the  past  —  e.g.,  Ashurst  and  Barr  (1983),  Ghoniem  et  al.  (1980,  1983),  Pope  (1990)  — 
that  has  significant  potential  to  treat  practical  turbulent  premixed  flames  due  to  its 
computational  efficiency. 

3.2  Experimental  Methods 

Apparatus.  The  fan-stirred  combustion  chamber  was  developed  by  Groff  (1987) 
based  on  an  original  concept  of  Semenov  (1965).  A  similar  arrangement  was  used  by 
Abdel-Gayed  et  al.  (1984).  The  chamber  is  quasi-spherical  with  a  260  mm  crossectional 
diameter  at  the  center.  Optical  access  was  provided  by  two  92  mm  diameter  windows  in  the 
end  wall  and  two  10  mm  diameter  windows  in  the  side  walls  (each  pair  of  windows 
mounted  opposite  one  another).  The  isotropic  turbulent  field  was  generated  by  four  fans 
mounted  at  90°  intervals  along  the  periphery  of  the  chamber  with  the  fan  flows  directed 
toward  the  walls.  Fansler  and  Groff  (1990)  show  that  this  arrangement  provides  an  istropic 
flow  field  in  the  central  region  of  the  chamber. 

The  combustible  mixture  was  prepared  by  filling  the  chamber  with  hydrogen, 
nitrogen  and  air  to  appropriate  partial  pressure  levels  and  then  mixing  the  gas  with  the 
fans.  The  mixture  was  ignited  by  a  spark  at  the  center  of  the  chamber.  Measurements  were 
limited  to  the  period  where  the  pressure  increase  within  the  chamber  was  small  so  that  the 


flames  corresponded  to  constant  pressure  deflagration  waves  in  an  unreacted  turbulent  gas 
having  constant  properties. 


Instrumentation.  Measurements  involved  laser  velocimetry  to  find  the  flow 
properties  of  the  unreacted  gas  and  flame  tomography  to  find  flame  surface  statistics.  Both 
measurements  involved  adding  small  kerosene  drops  (nominal  diameter  less  than  1  pm)  to 
the  unreacted  gas;  these  drops  disappeared  at  the  flame  surface  as  discussed  earlier. 

The  laser  velocimetry  system  was  a  dual-beam  forward-scatter  arrangement  that 
provided  both  single-and  two-point  velocity  statistics.  Both  measurements  used  frequency 
shifting  to  control  directional  bias  and  ambiguity  with  the  laser  beams  rotated  to  measure 
velocity  components  in  orthogonal  (vertical  and  horizontal)  directions.  The  single-point 
measurements  were  straight-forward  and  involved  directing  the  laser  beams  and  observing 
the  measuring  volume  (in  the  forward  scatter  direction)  through  the  large  windows.  The 
two-point  measurements  involved  directing  the  laser  beams  at  a  small  angle  through  the 
small  windows  to  provide  a  measuring  volume  that  was  70  mm  long.  Two  points  within 
this  measuring  volume  were  observed  with  two  traversable  detectors  through  the  large 
windows  to  provide  two-point  measurements. 

The  flame  tomography  measurements  were  similar  to  the  jet  flames.  A  pulsed-dye 
laser  was  used  to  form  a  thin  light  sheet  (200  pm  thick)  directed  through  one  of  the  small 
side  windows.  The  light  sheet  was  observed  normal  to  its  plane  from  the  one  of  the  large 
windows  using  a  35  mm  SLR  camera  whose  exposure  was  controlled  by  the  short  duration 
of  the  laser  pulse  (1-2  ps).  The  boundary  between  the  unbumed  and  burned  gas  could  be 
observed  due  to  the  absence  of  scattering  from  drops  within  the  burned  gas.  This  boundary 
was  analyzed  similar  to  the  turbulent  jet  flames  to  find  flame  surface  statistics. 

Test  Conditions.  Test  conditions  for  freely-propagating  flames  are  summarized  in 
Table  3.  All  tests  were  at  0  =  1.8,  which  is  the  neutral  preferential-diffusion  condition. 
Values  of  Sl  were  varied  by  varying  the  degree  of  dilution  with  nitrogen  to  yield  u’/Sl  in 
the  range  0.48-1.60.  The  turbulence  was  well  developed  with  turbulence  Reynolds 
numbers  in  the  range  1965-4195,  representative  of  high  intensity  practical  premixed 
turbulent  flames.  Characteristic  flame  thicknesses  were  generally  an  order  of  magnitude 
smaller  than  the  Kolmogorov  scale  of  the  turbulence,  which  is  characteristic  of  the  thin 
flamelet  regime. 

3.3  Theoretical  Methods 

Turbulent  flame  propagation  was  numerically  simulated  using  the  flame  advection 
and  propagation  algorithm  from  the  MIMOC  routine  (Ghoniem  et  al.,  1983),  combined 
with  a  statistical  time  series  simulation  of  flow  velocities  in  the  unbumed  gas.  As  noted 
earlier,  MIMOC  can  only  provide  a  two-dimensional  time  dependent  simulation  of  flame 
surface  motion.  This  is  a  major  limitation  but  it  was  desired  to  evaluate  the  simplified 
simulation  before  extending  it  to  treat  three-dimensional  effects. 

Other  major  assumptions  of  the  approximate  numerical  simulation  are  as  follows: 
constant  pressure  deflagration  wave;  infinitely-thin  flame  sheet  with  constant  unbumed  and 
burned  gas  scalar  properties;  stationary  homogeneous  isotropic  turbulence  in  the  unbumed 
gas,  unaffected  by  the  presence  of  the  flame;  neutrally-stable  flame  with  negligible  effects 
of  quenching,  i.e.,  relative  to  the  gas,  the  flame  propagates  normal  to  its  surface  at  a 
constant  laminar  burning  velocity;  and  mean  velocity  field  in  the  unbumed  gas  was  found 
assuming  that  the  flame  acts  like  a  spherically  symmetric  volumetric  source.  The  constant 
pressure  and  thin  flamelet  assumptions  are  conditions  of  the  experiments.  The  assumption 
that  the  turbulence  properties  of  the  unbumed  gas  are  not  affected  by  the  flame  is  an  open 


issue,  although  existing  measurements  suggest  that  the  volumetric  expansion  of  flame 
passively  convects  the  turbulent  field  of  the  unbumed  gas  away  from  the  ignition  source 
(Kwon  et  al.,  1991).  Neutral  stability  is  a  condition  of  the  experiments  while  effects  of 
quenching  are  not  large  because  the  present  flames  were  well  away  from  extinction 
conditions  and  had  modest  values  of  u’/Sl-  Finally,  the  volumetric  source  approximation 
yielded  an  excellent  representation  of  mean  velocities  in  the  unreacted  gas  measured  by 
laser  velocimetry  during  turbulent  flame  propagation  (Kwon  et  al.,  1991). 

Table  3.  Summary  of  Turbulent  Freely- Progagating  Flame  Tests3 


02/(02  +  N2) 

0.210 

0.210 

0.150 

0.150 

0.125 

u/SL 

0.48 

0.90 

0.80 

1.60 

1.00 

u’(m/s) 

1.2 

2.4 

1.2 

2.4 

1.2 

Sl  (m/s)b 

2.5 

2.5 

1.5 

1.5 

1.2 

Rej  =  u'o  Uuo/Vo 

1965 

3930 

2095 

4195 

2320 

Qk  (|im) 

21 

12 

13 

8 

12 

8l  =  Oq/Sl  (Jim) 

1.8 

1.8 

2.8 

2.8 

3.3 

3  Hydrogen,  air  and  nitrogen  mixtures  at  neutral  preferential-diffusion  conditions  (<>  =  1.8), 
in  the  thin  flamelet  regime,  at  3  atm.  and  298  ±  3K.  The  longitudinal  integral  length  scale  of 
the  turbulence  was  12.5  mm. 
bFrom  Lewis  and  von  Elbe  (1961). 

The  statistical  simulation  of  the  velocity  field  of  the  unbumed  gas  was  similar  to  the 
approach  discribed  for  the  jet  flames.  The  computations  simulated  mean  velocities, 
Gaussian  PDF's  of  velocity  fluctuations  and  temporal  and  spatial  correlations.  Cross 
correlations  are  zero  and  the  locally-homogeneous  approximation  is  not  needed  for 
isotropic  turbulence.  The  temporal  and  spatial  correlations  were  approximated  as 
exponential  functions:  this  is  reasonably  accurate  for  temporal  and  longitudinal  spatial 
correlations  but  does  not  represent  the  negative  portion  of  the  Frenkiel  function  shape  of  the 
transverse  correlation  (Hinze,  1975;  Tennekes  and  Lumley,  1972).  Nevertheless,  the 
approximation  was  adopted  for  the  transverse  spatial  correlation  as  well,  because  the 
significant  portion,  where  the  values  of  the  correlation  are  near  unity,  is  approximately 
exponential.  Similarly,  because  the  high  correlation  region  is  of  major  interest,  correlations 
were  decomposed  as  products  of  correlations  along  each  coordinate  direction  as  discussed 
in  connection  with  the  jet  flame  simulation. 

Under  these  approximations,  the  simulation  of  velocity  fluctuations  at  a  particular 
point  only  involves  the  seven  nearest  neighbors  of  the  point  in  space  and  time  for  a  two- 
dimensional  time-dependent  simulation  (the  15  nearest  neighbors  for  a  three-dimensional 
time-dependent  simulation).  This  arrangement  provides  a  Markov-like  process  that 
simulates  the  decay  of  correlations,  taken  to  be  exponential  functions,  for  all  space  and 
time.  Kwon  (1991)  discusses  evaluation  of  numerical  accuracy  of  the  method:  for  4000 
realizations,  r.m.s.  fluctuations  were  accurate  within  5%  and  temporal  and  spatial 


correlations  were  accurate  with  30%  (within  two  integral  scales).  A  larger  number  of 
realizations  would  reduce  the  uncertainties  of  simulated  properties. 


A  typical  realization  of  the  fluctuating  velocity  field  for  conditions  where  no  flame  is 
present  is  illustrated  in  Fig.  4.  The  effect  of  allowing  for  spatial  and  temporal  correlations  is 
evident  from  the  cell-like  structure  with  cell  size^  comparable  to  integral  scales;  in  contrast, 
simulation  of  only  probability  density  functions  yields  a  nonphysical  random  structure  with 
points  close  to  one  another  having  velocities  in  different  directions.  In  fact,  the  simulation 
of  Fig.  4  bears  a  striking  resemblance  to  the  instantaneous  turbulent  velocity  fields 
measured  by  Reuss  et  al.  (1990)  using  particle-image  velocimetry. 

3.4  Results  and  Discussion 

Unbumed  Gas  Properties.  The  turbulence  properties  of  the  unburned  gas  were 
measured  for  all  conditions  used  during  testing.  These  measurements  were  completed  in  the 
absence  of  flame  propagation.  Within  ±30mm  from  the  center  of  the  chamber,  the  results 
showed  mean  velocities  less  than  10%  of  velocity  fluctuations,  v'/u'  generally  within  10% 
of  unity,  and  variations  of  u'  and  v’  with  position  less  than  10%.  Velocity  fluctuations 
were  proportional  to  the  fan  speed  and  relatively  independent  of  the  pressure.  Spatial  and 
temporal  correlations  approximated  exponential  functions,  except  for  the  crosstream  spatial 
correlation  at  large  distances  as  discussed  earlier.  The  spatial  integral  scale  based  on 
streamwise  fluctuations  was  12.5  mm,  relatively  independent  of  position  and  fan  speed, 
within  ±  30  mm  from  the  center  of  the  chamber.  Temporal  integral  scales  were  inversely 
proportional  to  the  fan  speed. 

Measurements  of  mean  radial  velocities  in  the  flow  were  obtained  by  ensemble 
averaging  single-channel  measurements  made  in  the  presence  of  propagating  flames.  These 
results  yielded  a  excellent  correlation  with  estimates  for  potential  flow  from  a  volumetric 
source,  as  follows  (Kwon  et  al.,  1991): 

u  /  (ST  (pu  /  Pb  -  1 ))  =  (Tf /r)2  (1) 

Results  of  Eq.  (1)  as  well  as  the  direct  measurements  of  turbulence  properties  were  used  in 
the  simulation  to  estimate  flame  surface  properties. 

Flame  Surface  Statistics.  Typical  flame  tomographs  for  the  present  neutrally  stable 
conditions  at  u/Sl  =  1-6,  which  is  the  highest  value  considered,  are  illustrated  in  Fig.  5. 
The  tomographs  are  for  rf  =  15,  30  and  45  mm.  The  results  show  a  progressive  increase  of 
flame  surface  distortion  with  increasing  radius  (or  time  of  propagation).  This  behavior  is 
analogous  to  the  progressively  increasing  distortion  of  the  flame  surface  of  premixed 
turbulent  jet  flames  with  increasing  distance  from  the  flameholder  that  was  discussed 
earlier. 


Simulated  flame  surfaces  at  u'/Sl  =  1.6  (the  experimental  value)  and  3.2  also  are 
illustrated  in  Fig.  5.  Both  simulations  represent  the  trend  that  flame  surface  distortion 
increases  with  mean  radios  (time).  The  main  effect  of  increasing  u'/Sl  is  to  increase  the 
degree  of  distortion  of  the  surface  at  particular  values  of  rf  so  that  the  larger  u'/Sl  yields  a 
more  irregular  surface  with  Finer-grained  distortion.  Both  simulations  are  qualitatively 
similar  to  the  flame  tomographs,  however,  the  small-scale  distortions  are  better  represented 
by  the  results  for  u’/Sl  =  3.2  (which  is  twice  the  experimental  value).  This  is  reasonable, 
because  a  two-dimensional  simulation  is  expected  to  underestimate  effects  of  true  three- 
dimensional  turbulent  distortion.  More  quantitative  assessment  of  the  simulations  will  be 
undertaken  in  the  following. 


Fig.  4  T  v :  i  c  ,1 1  ^  l  mu  lotion  \  veloci 
i  n  i  t re  n  i  c  turhu  1  <?t~p . 


luct’.iriiicn 


Measured  and  simulated  values  of  the  fractal  dimension,  D3  =  D2  +  1 ,  are  plotted  a:; 
a  function  of  7 f  in  Fig.  6.  At  small  7f.  D3  is  nearly  2.0  which  is  representative  of  a  smooth 
spherical  spark  kernel.  The  subsequent  increase  is  due  to  progressive  deformation  of  the 
flame  surface  by  turbulence,  somewhat  analogous  to  the  behavior  of  D3  for  turbulent  jet 
flames  with  increasing  distance  from  the  flame  holder.  However,  maximum  values  of  D3, 
are  not  large,  generally  less  than  2.15.  This  is  caused  by  the  limited  propagation  time  of 
present  tests,  which  only  corresponds  to  0.4- 1.3  times  the  integral  time  scale  or  distances 
of  propagation  into  the  unreacted  gas  mixture  on  the  order  of  one  spatial  integral  scale  (the 
larger  displacements  of  7f  being  caused  by  effects  of  volumetric  expansion).  While  the  D3 
illustrated  in  Fig.  6  show  monotonic  grow-th  in  the  early  stages  of  flame  grow  th  that  were 
observed,  it  is  expected  that  D3  would  eventually  approach  values  in  the  range  2. 3-2.4  that 
are  characteristic  of  passive  isoclinic  surfaces  in  isotropic  turbulent  fields.  However, 
measurements  at  larger  values  of  7f,  or  smaller  spatial  and  temporal  integral  scales,  are 
needed  to  confirm  this  behavior.  The  results  do  show,  however,  that  the  rate  of  increase  of 
D3  with  distance  increases  as  u'/Sl  increases.  Thus,  correlations  of  D3  as  a  function  of 
u  /Sl  have  been  found  for  flames  of  particular  sizes  (Mantzaras  et  ah,  1987)but  clearly  do 
not  account  for  the  important  effect  of  distance  from  the  point  of  ignition. 

Simulated  values  of  D3  in  Fig  6  also  exhibit  a  progressive  increase  with  increasing 
flame  radius.  However,  use  of  the  experimental  value  of  u'/Sl  yields  simulated  values  of 
D3  that  are  significantly  smaller  than  the  measurements  for  each  value  of  7f.  However, 
doubling  u/Sl  for  the  simulation  yields  a  reasonably  good  estimate  of  the  variation  of  D3 
with  7f.  A  probable  reason  for  this  deficiency  is  that  out  of  plane  distortions  of  the  flame 
surface,  which  should  contribute  to  its  irregularity,  are  suppressed  by  two-dimensional 
simulations. 

Another  measure  of  the  distortion  of  flame  surfaces  by  turbulence  was  found  by 
measuring  the  perimeters  of  the  flame  surfaces.  For  present  conditions,  flame  wninkles  are 
not  large  and  it  is  reasonable  to  assume  that  surface  properties  are  isotropic  so  that  the 
wrinkled  surface  can  be  related  to  the  perimeter  as  follows  (Wu  et  ah,  1991a); 

AT/AL  =  2(Pr/PL)- 1  (2) 

Additionally,  for  neutral  preferential  diffusion  conditions  and  modest  stretch  rates,  Aj/Al 
is  a  measure  of  Sj/Sl  as  noted  earlier. 

Measured  and  simu,-ued  values  of  Pj/Pl  are  plotted  as  a  function  of  rf  in  Fig  7. 
The  measurements  show  tnat  Pj/Pl  is  near  unity  for  small  values  of  rf,  which  is 
representative  of  the  nearly  smooth  spherical  spark  kernel  in  this  region.  Similar  to  D3, 
however,  Pt/Pl  progressively  increases  with  increasing  rf,  with  the  rate  of  increase  being 
more  rapid  as  u’/Sl  increases.  Whether  Pj/Pl,  and  thus  Sj/Sl,  eventually  approaches  an 
asymptotic  value,  even  though  D3  reaches  an  asymptotic  value,  is  r  major  open  issue  for 
premixed  turbulent  flames  —  tests  involving  longer  progation  times  are  needed  to  help 
resolve  this  important  property.  Associating  Pj/Pl  with  Sj/Sl  indicates  that  it  increases 
with  rf,  with  u’/Sl  largely  controlling  the  rate  of  increase,  for  laboratory  scale 
measurements  like  tiiose  illustrated  in  Fig.  7.  Thus,  correlations  of  Sj/Sl  solely  as  a 
function  of  u’/Sl  are  only  appropriate  for  particular  values  of  rf,  as  recently  shown  by 
Trautwein  et  al.  (1990). 

Because  Dj  and  Pj/Pl  are  closely  related  flame  surface  features,  the  comparison 
between  measured  and  simulated  values  of  Pt/Pl  in  Fig.  7  is  similar  to  results  for  D3  in 
Fig.  6.  In  general,  the  simulations  yield  the  correct  trends  of  Pj/Pl  with  increasing  7f  and 
u’/Sl  hut  quantitative  agreement  is  only  achieved  by  using  values  of  u’/Sl  roughly  twice 


Observed  and  simulated  flame  surface  images  for  freely-propagating  flames: 
O2AN2+O2)  =  0.15;  r f  =15,  30  and  45  mm.  From  Kwon  et  al.  (1991). 


the  experimental  values.  As  before,  the  limitations  of  a  two-dimensional  simulation  are  the 
most  probable  cause  of  this  difficulty. 

4.  CONCLUSIONS 

A  theoretical  and  experimental  study  of  the  properties  of  turbulent  premixed 
hydrogen/oxygen/nitrogen  flames  was  completed,  considering  both  jet  and  freely- 
propagating  flames  at  turbulence  Reynolds  numbers  in  the  range  280-4200.  The  main 
canclusions  of  the  study  can  be  summarized  as  follows: 

1 .  Effects  of  preferential  diffusion  are  significant  for  high  Reynolds  number  flames, 
acting  to  reduce/increase  rates  distortion  of  the  flame  surface  (and  correspondingly  the 
local  turbulent  burning  velocity)  for  stable/unstable  conditions.  Behavior  at  unstable 
conditions  is  of  greatest  interest  for  hydrogen/air  or  oxygen  flames,  used  for 
propulsion  or  flame  cutting  applications,  because  lean  or  stoichiometric  mixture  ratios 
are  used  for  these  applications  which  is  within  the  unstable  regime.  In  contrast, 
behavior  at  stable  conditions  is  of  greatest  interest  for  heavy-hydrocarbon/air 
flames, used  for  propulsion  or  automotive  applications,  because  the  relevant  lean 
mixture  ratios  fall  in  the  stable  regime  for  these  fuels. 

2.  Spatial  variations  of  flame  surface  properties  were  quantified.  As  the  distance  from 
the  flameholder  or  point  of  ignition  increased,  even  though  unbumed  gas  turbulence 
properties  were  relatively  constant,  the  following  properties  increased  in  value:  flame 
surface  area  normalized  by  the  mean  flame  area  (and  thus  the  local  turbulent  burning 
velocity),  flame  brush  thickness,  and  the  fractal  dimension  of  the  flame  surface. 
However,  at  large  distances,  the  fractal  dimension  approached  values  appropriate  for 
isoclinic  surfaces  in  isotropic  turbulence,  e.g.,  D3  =  2. 3-2.4,  while  approach  to  the 
flame  tip  limits  the  maximum  value  of  the  other  properties  for  jet  flames.  These 
effects  imply  that  turbulent  premixed  flames  are  geometry  specific  while  the  spatial 
variations  imply  that  models  or  correlations  of  turbulent  flame  properties  based  on 
local  properties  alone,  which  don't  specifically  account  for  distance  from  the 
flameholder  or  point  of  ignition,  are  intrinsically  incomplete,  i.e.,  there  is  no  basis  for 
use  of  the  local  conditions  hypothesis  for  the  combustion  rates  of  turbulent  premixed 
flames. 

3.  Increasing  values  of  turbulent  fluctuations  relative  to  the  laminar  burning  velocity, 
u'/Sl,  tend  to  increase  the  rate  of  distortion  of  the  flame  surface  with  distance  from 
the  flameholder  or  point  of  ignition  and  also  increases  levels  of  fine-grained 
roughness  of  the  flame  surface.  Thus,  correlations  of  turbulent  burning  velocity  and 
flame  surface  fractal  dimension  solely  as  function  of  u'/Sl  can  be  found  because 
they  reflect  this  rate  of  increased  over  limited  ranges  of  flame  dimensions. 
Nevertheless,  such  correlations  are  incomplete  unless  they  account  for  development 
of  flame  surface  properties  with  distance  from  the  flameholder  or  point  of  ignition. 

4.  Estimates  of  turbulent  jet  flame  properties  based  on  a  contemporary  turbulence  model 
due  to  Cant  and  Bray  (1988),  which  allows  for  effects  of  flame  stretch,  yielded  only 
fair  agreement  with  measurements.  The  difficulties  are  attributed  to  effects  of 
preferential  diffusion  and  finite  laminar  burning  velocities  which  are  not  considered  in 
the  model,  and  the  use  of  the  local  conditions  hypothesis  which  is  intrinsically 
incomplete,  as  noted  earlier. 

5.  The  present  two-dimensional  numerical  simulation  successfully  predicted  the  trends 
cf  flame  properties  with  distance  from  the  flameholder,  or  point  of  ignition,  and 
u /Sl-  Quantitatively  accurate  predictions  near  the  base  of  the  jet  flames,  or  for  freely 


;u 


propagating  flames  having  diameters  larger  than  the  integral  length  scale,  could  be 
obtained  by  the  artifice  of  doubling  uVSl  in  order  to  compensate  for  the  limitations  of 
a  two-dimensional  simulation  of  three-dimensional  turbulence.  However,  the 
generality  of  this  artifice  is  unknown  and  it  is  less  effective  where  motion  of  the  flame 
as  a  whole  is  significant,  i.e.,  near  the  tip  of  jet  flames  or  when  the  flame  kernel  is 
small  for  freely  propagating  flames. 

6.  Three-dimensional  time-dependent  simulations  are  needed  to  definitively  evaluate  the 
effectiveness  of  the  present  numerical  simulations.  Fortunately,  extension  to  three 
dimensions  is  tractable,  even  for  practical  high  Reynolds  number  flames,  due  to  the 
computational  efficiency  of  the  present  velocity  simulation.  Issues  of  modification  of 
unhurried  gas  turbulence  properties  by  the  flame,  quenching  and  preferential  diffusion 
still  must  be  assessed.  The  present  stochastic  simulation  provides  a  promising 
framework  to  consider  these  effects  because  the  instantaneous  flame  surface  is 
simulated,  from  which  local  flame  curvature,  stretch  and  flow  acceleration  can  be 
found. 
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APPENDIX  B:  WU  ET  AL.  (1991a) 


of  Turbulent  Premixed  Hydrogen/Air  Flames  Turbulent  premixed  flames  have  numerous  practical  applications,  eg,  spark- 


40 


u 


kHz,  using  a  LeCroy  a/d  convener  (12  bit).  Power  spectra  were  found  from  64k  sample:  Kolmogorov  scales,  fit,  were  estimated  as  fiuo/(  u  oavgfiuo/vu) 1/4  from  Tennekes  and 

at  each  location  using  an  FFT  algorithm  and  an  lBM-AT/clone  computer.  Experimenta  Lumley  (1972):  they  are  generally  comparable  or  larger  than  the  flame  thickness  which  is 

uncertainties  (95%  confidence)  were  dominated  by  finite  sampling  times  and  are  estimate!  representative  of  the  thin  flamelet  regime, 

to  be  less  than  5%  for  mean  and  fluctuating  velocities  in  the  streamwise  direction  anc 
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Normalized  probability  density  functions  (PDF's)  of  streamwise  and  crosstrearr  50  dial  exponential  autocorrelations,  which  simplify  stochastic  simulations  of  the  flow  held 

velocity  functions  for  various  points  within  the  unreacted  gas  are  plotted  in  Figure  3  tBox  and  Jenkins,  1976),  would  be  a  reasonable  approximation  for  most  of  the  flow. 

These  results  are  for  the  same  flame  as  Figures  1  and  2,  however,  they  are  typical  of  othei  ,emPoral  integral  scales,  xu  and  tv,  are  plotted  in  Figure  6  for  position- 

conditions.  In  spite  of  the  anisotropy,  the  velocity  fluctuations  are  seen  to  approximate  corresponding  to  the  spectra  of  Figure  5.  The  streamwise  integral  lengtn  scale,  l!u,  is  also 
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yiim>  Position  Mean  flame  positions,  ff  and  the  r  m.s,  tluclnation  of  flanu  Flame  Perimeter  and  Area  A  measure  ol  the  distortion  ol  the  flame  surface 

position  or  flame  brush  thickness,  r'f.  are  plotted  as  a  function  of  height  above  the  burnci  was  found  by  measuring  the  perimeters  of  the  flame  surfaces  on  die  vertical  and  horuonial 
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hydrogcn/methane/air  flames  exhibited  faster  growth  rates  of  Sj/Sl  w‘th  increasing  jn  orcjcr  iq  compare  present  results  with  other  work,  where  the  assumption  of  isotropic 

u  o^vg/SL  ter  ncutraVunstable  conditions  than  for  stable  conditions  The  results  illustrated  surface  properties  is  generally  invoked  (Gouldm,  11>K7.  l.ovejoy,  1982.  Mant/aras  et  al  . 


smooth  geometrical  surface  (Mandelbrot,  1982).  This  follows  since  the  point  of  aitachmer  maximum  values  lhai  I):  can  achieve  Other  trends  wall  respect  lo  u',,jvg/S[.  alone  tend  to 

of  the  flame  ts  fixed  along  the  periphery  of  the  burner  passage  exit  and  the  turbulence  doe  be  obscured  by  limited  experimental  accuracy  and  effects  of  preferential  diffusion. 


surface  as  effects  of  turbulence  begin  to  dominate  the  combustion  process.  Effects  of  flamt  increases  with  increasing  distance  from  ihe  burner  exit  with  maximum  values  largely 

position  are  more  complex  due  to  the  boundary  conditions,  with  Xp/Huo  increasing  near  the  governed  by  the  length  of  the  flame  '1116  relative  turbulence  intensity,  u'lMvg/Si.  primarily 
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are  only  penmen!  for  similar  flame  conditions.  Again  limning  considerations  to  neutra  ra(e  of  mcreaic  0f  llame  disionton  for  stable/unstable  conditions,  similar  to  the  predictions 

conditions,  a  finite  relative  velocity  of  a  fluid  point  with  respect  to  'he  tlow,  like  SL  of  C|avjn  and  wilhams  ,1919  i9H2>  for  premixed  flames  in  weak  turbulence.  When 


Spatial  variations  of  flame  surface  properties  were  quantified  As  the  distance  from 


correspondingly  the  local  turbulent  burning  velocity)  for  stable/unstable  conditions. 
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Nomenclature  p  density 

Ai/At  average  surface  area  of  time-averaged  flame  surface,  of  wnnkled  interface  T‘  integral  lime  scale  based  on  velocity  component 
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APPENDIX  C:  KWON  ET  AL.  (1991) 


Flame  Surface  Properties  of  Premised  Flames  in  Isotropic  Turbulence:  G  value  of  g(Ax)  =  g(Ay) 

Measurements  and  Numerical  Simulations  Kolmogorov  length  scale 
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Introduction  flame  surface  development,  will  be  useful  for  developing  and  evaluating  models  of  the 

An  experimental  and  theoretical  investigation  of  free  turbulent  premised  flame  procesi. 

propagating  in  isotropic  turbulence  at  neutrally-stable  preferential  diffusion  conditions  i  Pas*  attemPB  *°  develoP  mode's  or  simulations  of  premixed  turbulent  flames  in  the 


Hydrogen  tod  nitrogen  <99.95%  parity)  were  supplied  from  commercial  cylinders  Elimt  IflUlfltfaaiU  A  pulsed  dye  laser,  providing  0.6  J  of  light  per  pulse  at 

while  dry  air  (dew  point  less  than  240  K)  was  obtained  from  laboratory  supplies.  The  5 14  5  11111  »  2  ps  pulse  duration,  was  used  for  the  FT  measurements.  The  laser  beam 

proper  partial  pressures  of  hydrogen,  nitrogen  and  air  were  mixed  together  with  the  fans  was  Wlt*1  spherical  and  cylindrical  lenses  and  directed  through  one  of  the  small 
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iithou,.i  drf/dt  at  *  particular  Ff  progressively  .ncreases  as  u'/Sl  increases  due  to  Simulated  values  of  Ff  exhibit  trends  similar  to  the  measurements  in  Figs.  8- 10, 

increased  levels  of  flame  suiface  distortion  by  the  turbulence  At  larger  umes  and  u'/Sl  namely,  rf  increases  with  lime  and  the  rate  of  increase  tends  to  increase  as  u'/Sl  becomes 

howeve  .  d7f /dt  increases  with  incr*  ong  time  as  well,  because  effects  ot  the  progressive  larger.  However,  the  simulations  generally  underestimate  Ff.  This  behavior  is  consistent 
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Table  2  Summary  of  Turbulent  Flame  Test  Conditions1 
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Fig  1  Orthogonal  velocity  fluctuations  m  the  unbumed  gas. 

Fig.  2  Temporal  power  spectral  densmes  of  velocity  fluctuations  in  the  unbumed  gas 
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APPENDIX  D:  WU  ET  AL.  (1991a) 
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Evaluation  of  Stochastic  Simulations  Kolmogorov  length  scale 

of  Premixed  Turbulent  Jet  Flames  name  length  based  on  d  mne  averaged  reaction  progress  variable  of  one-half 
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value  of  transverse  spatial  correlation  coefficient  at  Ax  avg  average  value 


urn  value  potential  to  provide  a  fundamental  description  of  all  aspects  of  premixed  turbulent  flames, 


ased  on  the  flame  propagation  algorithm  of  Ghoniem  et  al.  ( 14X3),  as  a  useful  first  st  streamwise  velocity  fluctuations  ami  integral  length  scales  at  the  center  of  the  burner  exit 

efore  committing  resources  to  extend  the  method  to  treat  three-dimensional  effects.  were  in  the  range  200- 1 100  faking  .he  characteristic  flame  ituckne  >s  to  be  6)  -  u„/S| 
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Data  from  four  flames  at  neutral  preferential  diffusion  conditions  were  used  t<  probability  density  (unctions  t 1*1)1*)  However,  velocity  fluctuations  *ere  anisotropic  with 

evaluate  the  simulations,  see  Table  I  for  a  summary  of  the  test  conditions  Marne  surface  u7  v’  in  ihc  range  1  2-3  0.  with  the  larger  values  observed  near  die  edge  of  the  flow 


similar,  except  for  the  appearance  of  a  local  maximum  at  low  frequencies  for  posmons  near  propagation  algorithm,  also  discussed  below,  requires  dial  each  point  on  ihe  (lame  surface 

the  axis  and  the  burner  exit,  Laufer  (1954)  observed  a  similar  spectral  feature  for  turbulent  ls  advccted  at  vel°cny  (u|t),  v(tn  and  propagates  normal  to  the  Home  surface  at  a  velocity 


8«J 


functions.  Cross  trc*m  temporal  correlations  approximate  exponential  (unctions  to  a  lesser  A  ,in,iJ  llullcr  involves  the  actual  specification  of  the  unhurried  gas  flow  piopertica 

degree  near  the  burner  exit,  however,  the  discrepancy  is  not  large  tn  the  region  of  the  (lame  'Hie  measurements  of  Wu  el  al  <  I  "9 1 )  cvere  limned  but  generally  indicated  dial  11, > a 

brush  (Wu.  1991).  properr't  -uriin  ihe  unhurned  gas  approximated  die  pi-  penies  ot  turbulent  pipe  thm 


op 
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Courant  condition  for  the  flame  propagauon  algorithm:  At  <  Ax  /  (2lulrnai),  see  Noh  and  twodimensional  simulation  implies  that  the  flame  is  i  ruled  surface  normal  to  the  plane  ot 

Woodward  (1976)  A  typical  run  time  for  4000  time  steps  on  a  DEC  5000  computer  was  ,he  lomo2I3Ph  which  reduces  the  potential  for  irrcgularu.es  of  surface  propcn.es  ... 

four  hours  comparison  us  die  true  three  dimensional  surf  ace 
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As  ar.hcipated  from  ihe  previous  results,  the  simulated  unreactedness  in  Figure  6  is  dimensions  should  be  tractable  even  tor  high  Reynolds  number  flames  due  to  die 

generally  greater  than  the  measurements  with  flame  lengths  overestimated  accordingly  computational  efficiency  of  the  velocity  simulation  I'hree  dimensional  lime 

Thus,  simulated  values  of  are  24  and  67%  too  large  for  u’o  aVj(/S|.  =  0  4  and  I  3,  and  dependent  results  are  particularly  needed  to  determine  whether  simulating  moments, 

u  J  u\  =  2,  the  latter  being  ratio  that  gives  the  best  simulation  of  flame  brush  thicknesses  PDF’s,  spatial  correlations  and  temporal  correlations  of  the  unburned  gas  is  adequate 


pressure:  d  -  11  mm.  =  3. 1  mm,  SL  =  3  5  m/s,  6L  =  17  pm. 
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